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ABSTRACT 

Aims. We report a detailed attempt to model the ortho-to-para abundance ratio of C-C 3 H 2 so as to reproduce observed values in the cores of the 
well-known source TMC-1. According to observations, the ortho-to-para ratios vary, within large uncertainties, from a low of near unity to a 
high of approximately three depending on the core. 

Methods. We used the osu.2003 network of gas-phase chemical reactions augmented hy reactions that specifically consider the formation, 
depletion, and interconversion of the ortho and para forms of the C-C 3 H 2 and its precursor ion C-C 3 H 3 . We investigated the sensitivity of the 
calculated ortho-to-para ratio for C-C3H2 to a large number of factors. 

Results. For the less evolved cores C, CP, and D, we had no difficulty reproducing the observed ortho-to-para ratios of 1-2. In order to 
reproduce observed ortho-to-para ratios of near three, observed for the evolved cores A and B, it was necessary to include rapid ion-catalyzed 
interconversion processes. 
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1. Introduction 

Cyclopropenylidene, C-C 3 H 2 , is a widely-distributed abun- 
dant organic ring molecule in the interstellar medium 

ijTTiaddeus et al. lll985[lMatthews & Irvinelli985UVrtilek et al.( 

Il987l) . Its “linear” isomer, propadienylidene (/-C 3 H 2 ), has 


lower abundance dCernicharo et al.M 19991 

Fosse et alJ 120011 

Tevssier et alJl2005l ICernicharo et al. Ill99ll iKawaeushi et alJ 

199lUTurner et alJl2000l). For each isomer. 

two equivalent H 


nuclei, of spin 1 / 2 , couple to generate ortho (nuclear spin = 
1 ) and para (nuclear spin = 0 ) species with spin statistical 
weights of 3 and 1, respectively. The nuclear spin wave func¬ 
tions of the ortho states are symmetric to exchange of the pro¬ 
tons whereas those of the para states are anti-symmetric. Since 
the Pauli exclusion principle requires total wave functions to be 
antisymmetric to exchange of protons, the symmetry of the ro¬ 
tational wave function of the molecule to exchange must be 
antisymmetric for ortho spin states and symmetric for para 
spin states if the ground electronic state is symmetric to ex- 
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change. The dependence of the exchange symmetry of the ro¬ 
tational state on the rotational quantum numbers depends on 
the structure of the molecule and on the symmetry of the elec¬ 
tronic state (iTownes & Schawlowlll975l) . Moreover, the con¬ 
version of molecular species from rotational levels with ortho 
spin functions to those with para spin functions does not occur 
efficiently if at all by radiative or non-reactive collisional mech¬ 
anisms in the gas, so that, in the absence of reaction (or strong 
binding with a surface), the sets of rotational levels with ortho 
and para spin functions can be regarded as distinct species. 


The cyclic isomer of C 3 H 2 has been observed in the in¬ 
terstellar medium in transitions belonging to both spin states: 
for 0 -C 3 H 2 , the observed lines include the b-typc transitions 
lio-loi, 2 i 2 - 1 oi, 3 i 2 - 3 o 3 , 32 i- 3 i 2 and 33 o- 32 i at 18.3, 85.3, 
83.0, 44.1, and 27.1 GHz respectively, whereas for / 7 -C 3 H 2 , 
the observed l 7 -type transitions include 2 o 2 -lii, 2 ii- 2 o 2 , 220 - 
2 11 and 37 2-3r3 lines at 82.1, 46.8, 21.6, and 84.7 GHz re¬ 
spectively llCer nicharO-^ '^Il999t Im atthews& lrvinelll98.5t 


Thaddeus et al. | |l'985[~|T ^akuwa^aLn20(HnMadd^iiet^ 

1989HMorisawa et alJ2005l) . With this many transitions, ortho- 
to-para ratios of reasonably high accuracy can be obtained, al- 
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lowing a detailed understanding of the chemical processes in¬ 
volved and possibly even the history of the source. 

In the terrestrial laboratory one often refers to two-types 
of mixtures for ortho and para species involving two protons: 
the common normal variety, in which the ortho-to-para ratio is 
three, and the equilibrium variety, in which the ortho-to-para 
ratio reflects thermal equilibrium. In the interstellar medium, 
however, the ortho-to-para ratio at low temperatures is proba¬ 
bly governed by kinetic rather than thermodynamical consider¬ 
ations, including processes on grains, and is likely to be time- 
dependent and may even reflect the original circumstances of 
the molecule formation 


ll.e Rourlolll2flflrih . In fact, detailed estimates of the steady- 
state interstellar ortho-t o-para H? ratio at 10 K show that it 
is approximately 10~^ jPagani et alJIlQQll iLe Bourl^ll99l[ 
Ike Bourlot et alJflQQOh . far higher than the equilibrium value 
of 3 X 10-^. 


Molecular hydrogen is unique in showing such a strong de¬ 
viation between thermal and normal values for the ojp ratio 
at temperatures in the vicinity of 10 K. For heavier species, 
this deviation is much smaller because more levels of each spin 
modification are thermally populated so that values of 3.0 
for the thermal ratio pertain to quite low temperat ures. For 
FI 2 CO, dPickens & Irvineiri 999t iKahane et al. Ill 98^ . this ra¬ 
tio is reached at temperatures as low as a; 15 K, while for 
C-C 3 H 2 , the ratio is reached under thermal conditions at tem¬ 
peratures under 5 K (iTakakuwa et alJEoOlh . Nevertheless, the 
actual ojp ratio in cold dense interstellar sources is determined 
by kinetic considerations and may well be significantly under 
three at 10 K clouds for heavy species with two equivalent pro- 
tons, as detected for the species FFCO, H 2 CS, and H 2 CCO 

dPickens & Irvinell999tlMinh et al.ll99ltl^ishi et alll99lll . 


Although the ortho-to-para abundance rat io of C-C 3 H 2 
had be en stud ied observationa ll y in T MC-1 bv iMadden et al.l 
dlMSl) and T^akuwa e t al.l d200ll) and in L1527 by 


ITakakuwa et alJd 200 lh. this work is motivat ed by the recent ob¬ 
servational work of iMorisawa et alJ d2005l) (hereafter, MFK), 
who measured the ratio in a more comprehensive manner for 
the six different cores within the TMC-1 ridge: A, B, C, P, 
CP, and E. Their measured ratios range from 1.4+0.7 (TMC- 
IP and CP) to 3.0+1.5 (TMC-IA and B). Note that Core P 
sometimes refers to some substructure near the CP condensa¬ 
tion, i n which case, it can be referred to as CP-b dHirahara et al.l 
Il992l) . but o ther authors have used Core P to be the same 
as Core CP dflartnuist et al. lEoOlh . which is the very well- 
studi ed core where the abundances of complex carbon chains 
peak dSmith et alj|2004l ). In additi on, TMCl-B is the so-called 
ammonia peak dHowe etaill 19961) . The current understanding 
is that some of the cores (e.g. Core B) are more evolved than 
others (e.g. Core CP) althoug h the cause of the differential evo¬ 
lution is not fully understood dFfartaiiist et al nmi owe et ID 
ll996HMarkwick et al. l 20 ^. 


Table^shows the o bserved 0 /p abundance ratios of cyclic 
C-C 3 H 2 in the six cores dMFKii along with measured fractional 
abundances (X) of 0 -C 3 H 2 with respect to H 2 and estimated 
ages of the cores. Although the determination of cloud ages is 
still contentious, there is a sufficient consensus concerning the 
chemical indicators of age to use the results in this paper. It can 


be seen that, as noted bv iMFKl the ortho-to-para ratio seems to 
be generally correlated with the evolution of the core; i.e., for 
all cores other than Core E, the older the core chemically, the 
higher the ortho-to-para ratio within large uncertainties, espe- 
cia lly for the higher ortho-to-para ratios. 

ImEKI used a simple model of a small number of gas-phase 
chemical reactions to understand their measured abundance 
ratios in terms of the evolution of the cores. They showed 
that the ortho-para chemistry of C-C 3 H 2 is rather simple 
compared with other heavy molecules that have been studied, 
so that a complete understanding of the ojp abundance ratio 
is possibly obtainable. Our detailed model is based on their 
simple one, with, as discussed below, the important reactions 
for the formation and depletion of the ortho and para forms of 
C-C 3 H 2 and its protonated ion precursor embedded in a large 
chemical network. 

In this paper, we report the results of our detailed gas- 
phase chemical model of the ortho-to-para abundance ratio of 
C-C 3 H 2 . The remainder of this paper is divided as follows. In 
Section 2 we discuss the details of our gas-phase chemical 
model, while in Section 3, we illustrate our results for the ortho- 
to-para ratio of C-C 3 H 2 as a f unction of time. We also compare 
our results with the model of lMEKl As will be seen, our stan¬ 
dard results are not capable of producing ojp ratios greater than 
2 . 1 , so we consider some additional direct chemical processes 
that can convert the para form of the ring molecule to its ortho 
form. With these processes, we are able to account for the pos¬ 
sibly high ojp ratio in the more evolved cores. We then subject 
our results to an uncertainty analysis. Section 4, the discussion, 
concludes the paper. 

2. The model 

The chemical reaction network used for this wo rk is an up¬ 
dated version of the osu.2003 gas-phase network jSmith et alJ 
12004 . Several versions of this network are now available at 
the URL http://www.physics.ohio-state.edu/~eric 
/research. html. Eor this work, we added several other mod¬ 
ifications. Eirst, and most importantly, we added reactions to 
account for the chemistry of both forms of C-C 3 H 2 and its pro¬ 
tonated precursor C-C 3 H 2 (see below). Note that this ion has 
three-fold symmetry so that its protons are all indistinguish¬ 
able: the ortho f orm has a t otal nuclear spin of 3/2, while the 
two para forms ilOkall2nn4 each have a total nuclear spin of 
1 / 2 , although only one of these exists for a given rotational 
state dTownes & Schawlowll975l) . Secondly, we added new ex¬ 
perimental information on the neutral products and total rate 
coefficients for dissociative recombination reactions involving 


several hydrocarbon ions dLarson et alJIlQQfl 

Semaniak et al. 

1998 

iKalhori et al.ll 2002 Hnerkatch et al.lll999 

iMitchell et al. 

2003 

). Thirdly, we ignored the existence of the linear isomer of 


lower abundance - /-C 3 H 2 - to simplify the problem for C-C 3 H 2 . 


2.1. Formation and destruction of C-C 3 H 2 

The formation of C-C3H2 begins with the radiative associa¬ 
tion of C3H^ with H2 to produce the precursor ion C-C3H3. 
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Table 1. Some C-C3H2 observations of the cores in the TMC -1 ridge 


Source 

Observed ojp Ratio'' 

Observed X(o-C 3 H 2 ) x 10 '' Estimated Age (yr) 

Cores with high ojp ratio ■ ■ ■ 




TMC-1 A 

-1.3 

7 5+^^ 

~ 10^ * 

TMC-IB 

3.0![;];3.0±0.16'' 

'•^-5.4 

~ lONlO^ * 

TMC-IE 

2.8!?:‘ 

^•^-3.4 

~ lO'^-lO^ 

Cores with low ojp ratio ■ ■ ■ 




TMC-IC 

1.8!°;’;2.4±0.ff 

6.0!*;^; 13 ' 

~ 10^ * 

TMC-1 CP 

1.4!°;«;2.56±0.31‘' 

57 

~ 10^ 6 X 10“* * 

TMC-1 D“ 

^ -^-0.7 


6 X 10'* * 


° [Hirahara et alHl992h . Core D is Core CP-b according to this source. 

'’ ImfkI unless noted. The fractional abundances X of o- C 3 H 2 are with respect to H 2 . Based on the MFK logarithmic observables, we 
carried out slightly d ifferent error calculat ions, because here A logX is not much less than 1 , so that two Icr random errors AX_ 36 AX+ 
should be calculatedj Wakelam et all200.'5h . 

‘‘ Maddetjil99(ir 

Cox et all 


f Saito et al. 1 

|2002 

) 

) 

^ Smith et alJ 

l2W4 

'' Hartauist et al. 1 J200lh 


This association is not the only reaction these two reactants 
undergo; the H-atom transfer channel to form c-CsH, + H 
may be competitive, alth ough it is unclear whether or not the 
reacti on is endothermicjjHerbstet^^ 1984 t iMaluendes et alJ 
ll 993 h.Tnthe OS U.200 3 jSmith et a 112004 and rate99 networks 
jLe Teuff et aklEoP^ the H-atom transfer is assumed to be a 
minor channel at 10 K. Once the cyclic C3H3 ion is formed, its 
major destruction is by dissociative recombination, one chan¬ 
nel of which leads to C-C3H2: 


C3H+ H- H2 ^ 

C-C3H3 + hv (major) 

(1) 


C-C3H2 -1- H (minor) 

(2) 

c-C3H^ -h e ^ 

C-C3H2 H- H 

( 3 ) 


c-CsH-h ... 

( 4 ) 


where the ... in reaction 0 refer to the fact that the additional 
products H -H H or H2 have not been explicitly determined by 
experiment. According to the recent experiment performed by 
lAngelova et al. I 1 I 2 OO 4 I) on the dissociative recombination of 
C3H3, product channels containing molecules with three car¬ 
bon atoms account for 90 . 7 % of the products. The experiment 
was not able to differentiate among the C3Hm species, nor to de¬ 
termine if linear or cyclic species were present, so the branch¬ 
ing ratio for the important products C-C3H2 -1- H remains a free 
parameter. 

Once produced, C-C3H2 is destroyed by reactions with a va¬ 
riety of ions while its precursor ion can be destroyed by neutral 


products CsH + ... products 



Fig. 1 . Dominant reaction diagra m for the standard models. 
Some processes not considered bv iMFKl are indicated by bold¬ 
face. 


species in addition to dissociative recombination: 


C-C3H2 -1- [HJ, HCO+, H3O+] 

^ c-C3Hj^ -H X 

(5) 

C-C3H2 -H [C+, S+, Si+] 

^ products 

(6) 

C-C3H+ -H [C, 0 , S, Si] 

^ products 

(7) 


Note that from now on, the cyclic designation c- will be under¬ 
stood and not written explicitly. 

One can see that in the absence of C3H2 reactions involving 
ions other than protonating ones, of ion-atom reactions involv¬ 
ing C3H3, and of dissociative recombination of this ion to form 
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any products other than C 3 H 2 , there is a quasi-cycle or loop 
in which once the protonated ion C 3 H 3 is formed, dissocia¬ 
tive recombination and protonation serve only to convert these 
species into one another. As these processes occur, ortho-para 
conversion also occurs, as we shall now discuss. 

2.2. Para-to-ortho conversion: standard and 
alternative models 

The dominant set of reactions involving the ortho/para chem- 
istry o f C 3 H 2 in our standard models, based on the model of 
ImfkI is listed in Table |2l along with additional reactions for 
two alternative models. The important reactions in the stan¬ 
dard approach are shown in Fig. ^ where the rate coefficients 
and their numerical subscripts are defined in the table. In Table 
| 2 j the reactions are on the far left, with the rate coefficients and 
branching fractions next to them. The branching fractions are 
written in terms of the rate coefficients, which are numbered 
accordin g to the ir position in the table. Reactions in the simple 
model of] ImfkI are indicated by a dark circle. 

The basic structure of the reaction scheme commences with 
the production of by radiative association because the 

H 2 reactant, as discussed in the Introduction, is overwhelm¬ 
ingly para in nature for cold clouds. Thus the zero nuclear 
spin of H 2 and the spin of the individual proton in C 3 H^ can 
combine only to form a resultant nuclear spin of 1 / 2 , which 
corresponds to a para form of C 3 H 2 . Once P-C 3 H 2 is formed, 
it undergoes dissociative recombination to form both p- and o- 
forms of C 3 H 2 (as well as other neutral products). The C 3 H 2 
species can then undergo protonation reactions with a variety 
of ions designated HX^; these processes lead to both ortho and 
para versions of the protonated ion. As the loop reactions con¬ 
tinue, sizable proportions of 0 -C 3 H 2 are built up. In competi¬ 
tion with the production of ortho species within the loop, other 
reactions serve to force the system to exit the loop, including 
(i) non-protonating ion-molecule reactions with C 3 H 2 , (ii) dis¬ 
sociative recombination of C 3 H 3 to form products other than 
C 3 H 2 , (iii) ion-atom reactions involving C 3 H 3 . Of course, in 
a complex model, some of the products formed outside of the 
loop can eventually lead to reactions such that the products are 
once again C 3 H 2 and its protonated ion. In addition to the cy¬ 
cle of protonation and dissociative recombination, there may 
also be direct conversion processes; such processes will be dis¬ 
cussed later but can be seen in Table|2las part of our alternative 
models. 

In order to determine the selection rules and branching ra¬ 
tios for ortho/para products in protonation and dissociative re¬ 
combination, we have used the approach of loO il20n4 . in 
which angular momentum algebra is applied to the addition of 
the nuclear spins of the reactants to determine the spin angular 
momenta of the intermediate complex, which then dissociates 
into product states of given spin angular momentum. In the lan¬ 
guage of the rotation group, if one considers the two reactants 
to have proton nuclear spins of I\ and I 2 , the combined spin 
state can be written as Sip ® Sp. This state leads to the in¬ 
termediate state Sp+p © ^ 7 ,+ 7,-1 © ■ ■ • © S\p^p_\, which means 
that the complex can have total proton spin angular moment 


ranging from I\ -H I 2 down to |/i - / 2 I. Each of these angular 
momenta can then dissoci ate into pro duct spin states given by 
angular momentum rules ilokJ2nn4 . The approach of Oka is 
most valid for processes in which many rotational states of the 
products can be produced. These processes can be exothermic 
reactions or even thermoneutral reactions if the energy is suf¬ 
ficiently high. In this limit, one need not worry much about 
the rotational quantum number J, since many different J states 
are produced. For low-temperature thermoneutral processes, on 
the other hand, where a small number of J states may be pro¬ 
duced, the role of the individual J states needs to be taken 
into accou nt in a compr e hensive quasi-phase-space theoretical 
appro ach 


122li. Gerlich & Schlemmerl 


120021 ) . Although the approach of lOkal J2004 lends itself natu¬ 
rally to the idea of a reaction intermediate in which all protons 
are indistinguishable from one another and can be facilely ex¬ 
changed during the life of the complex, it is also possible to 
remove this assumption. In fact, in our use of the treatment 
for the ion , we assume that the protonation reactions occur 
by simple transfer of a proton rather than formation of a com¬ 
plex consisting of C 3 H 2 a nd H/. Our branching fractions are 

who derived theirs with the 


the same as those used by llVih Kt who derived theirs will 
approa ch of lOiiackl ill 9771) . which is equivaleht to that of| 

il2nn4 . 


Ih our staudard models, we do hot cohsider direct (ioh- 
catalyzed) ihtercohversioh processes betweeh the ortho ahd 
para forms of C3H2 because competitive reactious with abuu- 
daut iohs that might catalyze such a couversiou also occur. For 
example, the reactiou of C3H2 with protous is likely to produce 
exothermic products such as 03 !!^ -H H 2 rather thau uudergo 
the exchauge of protous heeded to cause ortho/para chauge 
while the reactious with protouatiug ious such as HCO^ aud 
are likely to uudergo protouatiou. Nevertheless, we hud, as will 
be showu Sectiou lT^ that iu the abseuce of such iutercouver- 
siou/exchauge reactious, the calculated ortho-to-para ratio for 
C3H2 cauhot be large euough to explaiu observed values hear 
3.0. So, ih our alteruative models we also cousider iutercou- 
versioh of o- aud P-C 3 H 2 by iou-catalyzed exchauge processes 
with the four ious H^, H 3 , H 30 ^ aud HCO^, the latter three oc- 
casioually being designated by HX^. For these processes, the 
approach of lOkal ll2004 leads to the detailed balance expres¬ 
sion in the limit that kT exceeds the ortho-para energy differ¬ 
ence; 


k{p o) ^ ^ ^ 3 

k{o p) ki 


where i is an index of the four pairs of interconversion pro¬ 
cesses considered. We simplify the forward and backward pro¬ 
cesses such that the reactants produce a complex with a given 
rate coefficient k - ki + k^p the complex then dissociates into 
0 -C 3 H 2 on 3/4 of the collisions and / 7 -C 3 H 2 on 1/4 of the 
colli s ions r egardless of the initial reactants. In the theory of 
lOkal (l2004 . the relevant formulae are k{o p) - 116 k and 
k(p —> o) - 1/2 k for isomerization reactions involving a pro¬ 
ton. For the unknown rate coefficients ki+k^i, we choose a stan¬ 
dard value equal to the ion-dipole collision rate. To the best of 
our knowledge, there is no evidence for or against the supposi¬ 
tion that these particular ortho/para exchange reactions are not 
competitive with reaction. Perhaps the strongest evidence that 
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Table 3. Default initial conditions 


Low-metal abundance 

Other parameters 

He 

6 . 00 (- 2 ) 


N 

2.14(-5) 


0 

1.76(-4) 


Hz 

5.00(-l) 

Temperature = 10 K 

C+ 

7.30(-5) 

rifi = 2 x 10 '* cm“'* 

S+ 

8 . 00 (- 8 ) 

/-= 1.3x10^*^ s-‘ 

Si+ 

8.00(-9) 

C/O ratio = 0.41 

Fe+ 

3.00(-9) 

Ay = 10 mag 

Na+ 

2.00(-9) 


Mg+ 

7.00(-9) 

0.5 k 2 

C3H+ + e-4 C3H2 + H 

P+ 

3.00(-9) 

0.5 k 2 

- 4 C 3 H+... 

CP 

4.00(-9) 


e 

7.3107(-5) 



Note : a(-b) = a x 10 with respect to H abundance 
Table 4. Summary of variations of standard model parameters 


Parameters Varied 

Range 

1) C/O ratio 

0.21 - 1.20 

2) Metal depletion 

low-metal - additional deple¬ 
tion of real metals by factors 
of 2 

low-metal - high-metal 

3) Cosmic-ray ionization rate 

1.3x10-** S-* - 1.3x10-*® S-* 

4) Gas-density hh 

2 x 10 ** cm-* - 2 x 10 * cm-* 

5) C 3 H 2 + H fraction 

0.50- 1.00 


such processes occur is discussed in ICordonnieret m toodi . 
where the exchange mechanism via the complex Hj is mea¬ 
sured to occur in the reaction between H 2 and in competi¬ 
tion with a simple proton hop. 

2.3. Parameters and models: default & sets of models 

1-6 

In our calculations of the ojp abundance ratio of C 3 H 2 , we start 
with a default set of initial conditions, listed in Table 0 used 
with our standard model of reactions. The default model con¬ 
tains so-called low-metal abundances (depleted heavy elements 
with C/O = 0.41) with the additional assumption that the initial 
abundances are atomic except for hydrogen, which is molecu¬ 
lar. In addition, the total hydrogen density is set at uh = 2 x 10"^ 
cm“^, the temperature at 10 K, the visual extinction at 10 mag, 
and the cosmic-ray ionization rate ^ at 1.3 x 10^'^ s“*. Finally 
the unknown branching fraction for the dissociative recombi¬ 
nation of C 3 H 3 is set so that the products C 3 H 2 -H H and C 3 H 
-I- ... are produced in equal amounts. The chemistry is followed 
with a single-point model under constant physical conditions. 

In an attempt to determine how the computed ofp abun¬ 
dance ratio for C 3 H 2 vs time depends upon an assortment of 


parameters, we have varied the C/O elemental ratio, the de¬ 
gree of metallicity, the cosmic ray ionization rate, and the gas 
density, as well as the neutral product branching ratio of the 
dissociative recombination reaction of the ion C-C 3 H 3 . These 
parameters affect how long the C 3 H 2 ortho-para system can 
stay in the loop shown in Fig. ^ Table 0] shows the range of 
the parameters considered in this work for our standard list of 
additional reactions. For each parameter, a number of different 
runs was attempted within the ranges listed in the table. The 
variation of parameters with our alternative model is more re¬ 
strictive, as will be discussed in Section lT^ 

3. Results 

3.1. Default and other standard sets of models vs 
observation 

The calculated ojp ratio for C 3 H 2 as a function of time in our 
default model is shown in panel a) of Fig.|2 Superimposed are 
the observational results for the cores in TMC-1; e.g.. Cores 
CP(D) and C, for which very short lifetimes are believed to 
pertain. Core A, with an intermediate lifetime, and Core B, for 
which a long lifetime of » 10^ yr is normally suggested. Core 
E is included although its lifetime is rather uncertain; the high 
ojp ratio su ggests some evolution , but the measured NH 3 /CCS 
ratio is low llHirahara et al II 1992 I) . suggesting a young chemi¬ 
cal age dMF^ Uncertainties in the lifetimes are large and not 
included in the figure. 

With our default model, we obtain the result that the ojp 
ratio is near unity at times relevant to Cores CP(D) and C, in 
reasonable if not excellent agreement with observed ratios, but 
never increases very much with time, reaching a steady-state 
value under 1.5, in some disagreement with the measured ra¬ 
tios of S.OC/j'j) and 3.0(/;[ j) for Cores A and B, where we have 
used our recalculated uncertainties (see footnote in Table 0. 
The calculated fractional abundance for 0 -C 3 H 2 is shown in 
panel c) of Fig.|2l where it can be seen that the agreement with 
the lower age cores is fair, with the model result being some¬ 
what high, is excellent for Core A, but is particularly poor for 
Core B. The calculated result is a typical one for pseudo-time- 
dependent models, in which most organic species have high 
abundance only at so-called “early time.” If the age of Core 
B were 10® yr, then the disagreement would be eliminated for 
the fractional abundance of 0 -C 3 H 2 although the disagreement 
with the observed ofp ratio would remain. Using the abund ance 
ratio of NH 3 to CCS as the criterion of chemical age, ImfkI con¬ 
sidered the age of core B to be comparable to Core A, but this 
analysis is unclear because the abundance ratio is indeed larger 
for Core B, although not by a huge factor. 

Given the probable disparity between the results of the de¬ 
fault model and observations for the ojp ratio at long times, 
we varied the five parameters discussed in Section 2.3 over the 
ranges listed in Table to obtain a variety of results with the 
standard sets of models 1-6. Note tha t the hi ghest density used 
is similar to the densities derived bv iMFKl: our default value 
is the more traditional one. The result was little change in the 
quantities of interest except when we varied the branching ra¬ 
tio for the products of the dissociative recombination branching 
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Fig. 2 . Calculated values of the ortho-to-para C3H2 abundance ratio (panels a and b) and fractional abundance of 0-C3H2 (panels 
c and d ) with resp ect to H? vs time compared wit h observed values in the cores of TMC- 1 . Note that the superimposed solid 
iIMFKI) and dotted llMaddenll 99 (ilCox et al.li^ error bars are from different observations. See text for explanation of models 
used. 


fraction. When the branching fraction for the products C3H2 + 
H is increased to unity, the calculated ojp ratio for C3H2 in¬ 
creases to about 1.9 at steady state. An increase in the ratio to 
about 2.1 can be obtained with the following other parameters: 
C/O elemental ratio = 0 . 2 , metal depletion a factor of two, no 
change in the value of and a higher density of 2 x 10^ cm“^. 
The results of this model, labeled our “best” standard model, 
for the ojp ratio and abundance of 0-C3H2 are shown in pan¬ 
els a) and c) of Fig. |2 The results for the ojp ratio are now 
barely within the large observational error bars, although the 
results for 0-C3H2 are still poor for the evolved core B, and 
would be poor even if the age of this core were reduced by 
an order of magnitude. As is well-known, the only way to in¬ 
crease the calculated abundance of carbon-chain type species 
at late times is to use a carbon-rich elemental abundance. With 
our best standard model changed only so that the elemental 
carbon-to-oxygen ratio C/O = 1 . 2 , we get the results shown in 
panels b) and d) of Fig. | 2 ] where it can be seen that the ob¬ 
served fractional abundance of 0-C3H2 as a function of core 
age is reasonably reproduced albeit somewhat high a t all times 
with respect to the most recent measurements liMFKh . although 
the 0 /p ratio has lost all of its time dependence. Nevertheless, 
the model results for the ojp ratio lie (barely) within the obser¬ 
vational errors except for Core B. 


3.2. Models with direct interconversion 

In order to reach an ojp ratio of 3.0 at steady-state, we 
included direct interconversion between the para and ortho 
forms of C3H2, as shown in Table |2 For the alternative ap¬ 
proach known as “Trial 2 ”, the conversion reactions are actu¬ 
ally assumed to be more rapid than protonation by a factor of 
five for the protonating ions HX^. The results of Trial 2 using 
the same set of parameters as used for the best standard model 
are also shown in panels a) and c) of Fig .|2 Here the calculated 
ratio reaches a steady-state value of near 2.7 and is in good 
agreement with observation over the complete time range. On 
the other hand, the agreement with the fractional abundance of 
0-C3H2 for late times (Core B) is once again poor, even if the 
age of this core is reduced to 10 ® yr. If, to remedy this latter 
deficiency, the C/O elemental ratio is set to 1 . 2 , then the frac¬ 
tional abundance of 0-C3H2 at late times is in fair agreement 
with the observed value in Core B but the calculated ojp ratio 
never rises above 2 . 1 . If the less dramatic assumption is made 
that the rate coefficients for interconversion are the same as for 
protonation (“Trial 1 ”), then the steady-state ojp ratio reaches 
a value of 2 . 5 , also in reasonable agreement with observation 
in the older cores. But the same problems occur with the frac¬ 
tional abundance of 0-C3H2 as occur for Trial 2 . Models with 
a C/O elemental ratio of 1.0 are no more successful. We con¬ 
clude that no simple pseudo-time-dependent model of the type 
considered here explains all of the data. 
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3.3. Comparison with simpie modei 

ImfkI developed a simple model to calculate the 0//2 C 3 H 2 
abundance ratio. In their model only a few types of reactions 
are considered (see • marks in Table Ell. In the limit of steady 
state, the ratio can be shown to be 


nio-CjWi) 

n(p-CiW2) 


(5 + 2,4>) X + 30 + 3 

(I +0 + 3 


3(0 + 1) 
0 + 3 


(9) 

where 0 is koRnie)Ikion-yniY), the ratio of the product of the 
dissociative recombination rate coefficient of € 311 : 5 ^, and 
the electron abundance, n{e), divided by the product of the rate 
coefficient for ion-atom reactions of C 3 H 3 , kion-y, and the abun¬ 
dance of reactive atoms, n{Y). The parameter 0 is essentially 
related to the ability of the reactants to stay in the loop shown 
in Fig.n If one makes the assumption that H 2 is overwhelm¬ 
ingly para in nature, then the ojp ratio of C 3 H 2 is determined 
solely by the parameter 0 , with a large value of 0 leading to an 
ojp ratio of 3 and a small value of 0 leadi ng to a ratio of 1. 
With the values for n{Y) and n{e) assumed bv lMFKl the former 
ratio pertains at steady state; their assumptions contain a rather 
high fractional abundance for electrons, which we would asso¬ 
ciate only with the high-metal case of elemental abundances. 
With the more standard low-metal case, the computed ratio be¬ 
comes very close to unity, in agreement with our more d etailed 
results. With the high-metal-type parameters chosen bv ImfkI 
their time-dependent calculation shows that the ojp ratio of 
C 3 H 2 shifts dramatically from 1 to 3 at some intermediate time. 

But there are two omissions in the theory leading to the 
simple steady-state formula for the ojp ratio: (i) destructive 
reactions involving C 3 H 2 with a host of non-protonating ions 
such as that make the system leave the loop, and (ii) dis¬ 
sociative recombination leading to products other than C 3 H 2 + 
H, which again have the effect of forcing the system to leave 
the loop. The effect of including both of these processes in the 
high-metal case makes the simple steady-state for the ojp ratio 
a more complex one and reduces the calculated ratio to val¬ 
ues close to those calculated in our more detailed approach. In 
other words, it is still difficult to produce large ojp abundance 
ratios in the absence of direct inter-conversion reactions. 


3.4. The roie of uncertainties in rate coefficients 

Another important aspect that we considered in our study is 
the uncert ainty in the rate coeffi cients. Using the method de¬ 
scribed in I Wakelam et al ] il2nn.5l) . we computed errors both in 
the fractional abundance of 0 -C 3 H 2 and in the 0 // 7 -C 3 H 2 abun¬ 
dance ratio for the default model. Briefly, this method consists 
of randomly generating new sets of reactions of the osu.2003 
database within their uncertainty range using a Monte-Carlo 
procedure. We ran 2000 different runs which gave 2000 values 
of the 0 -C 3 H 2 abundance and of the 0 /P-C 3 H 2 ratio. The error 
in the abundance and ratio is defined at each time as the enve¬ 
lope containing 66 % of the computed values (Icr). We chose 
a Icr error because the observations are also given an observa¬ 
tional error of Icr. We utilized an uncertainty of a factor of two 
for the rate coefficients of the reactions listed in TableEl For the 
reactions of this table that have more than one product channel, 
we kept the branching ratio constant, because they were derived 
from theory and, as a consequence, have no statistical error. 

We found a maximum error of +0.4 for the 0 // 2 -C 3 H 2 ratio 
at 10^ yr. This error decreases to +0.2 at 10^ yr before increas¬ 
ing again to +0.3 after 10® yr. As shown in Fig.|3] the computed 
uncertainty in the ratio is not large enough to significantly af¬ 
fect any of our conclusions. The computed error in the frac¬ 
tional abundance of 0 -C 3 H 2 increases with time. Specifically, 
the error domain is [Xo-CjHj/I-S, l.SXo-Cj/rJ at 10 "^ yr and 
[X 0 -C 3 H 2 / 4 , 4 X„-c 3 h,] at 10^ yr. In this case again (see Fig.|3}, 
the computed uncertainty does not improve the agreement with 
observation significantly. 


4. Discussion and Conciusion 

New observations bv iMFKi of the ortho-to-para abundance ratio 
for the cyclic isomer of C 3 H 2 in the six cores of TMC-1 indi¬ 
cate that this ratio is under 2.0 for two of three cores thought to 
be chemically young and is closer to the thermal value of 3.0 
for those cores thought to be chemically more evolved. The 
parameter used to determine chemical age is the abundance 
ratio of NH 3 to that of CCS. The ojp ratio is determined by 
a chemistry in which the cyclic ion C 3 H 3 is formed initially 
solely in its para state, but by a process of dissociative recom¬ 
bination and re-protonation, a significant ortho/para abundance 
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ratio, near the thermal value of 3.0, can be built up given suf¬ 
ficient time for high-metal abundances if comp etitive reactions 
are not dominant. Although a simple model bv iMFKl indicates 
agreement with their observations for five out of six cores if 
high-metal abundances are used, our more detailed calculations 
with an extended version of the osu.2003 network of reactions 
do not give such a definitive result, mainly because of compet¬ 
itive processes. In fact, with standard parameters varied over 
rather wide ranges, we find it difficult to obtain o/p ratios ex¬ 
ceeding 2.0 at any time. 

We are able to calculate higher o/p ratios for C 3 H 2 only if 
we assume that the ortho and para forms of this species can 
be converted into one another by exchange processes involv¬ 
ing and protonating ions HX^. This assumption is based 
on some experiments involving bv ICordonnier et ahl 1 I 2 OOOI) 
but here runs up against the problem that both and HX^ will 
likely react with C 3 H 2 and the competition between reaction 
and ortho!para conversion may well favor the former. Thus, 
there is still no definitive chemical proof that gas-phase pro¬ 
cesses can produce ortho jpara ratios for cyclic C 3 H 2 near the 
thermal value of 3.0. If, then, the observational values in this 
range are sufficiently accurate to distinguish correctly between 
low and high ratios, either the exchange reactions we consider 
between 0 -C 3 H 2 and P-C 3 H 2 do occur and are actually favored 
over reaction, our branching fractions or other aspects of the 
gas-phase chemistry are in error, or one must consider gas- 
grain interactions. Future observations with a lower uncertainty 
will help to determine which if any of these choices is the cor¬ 
rect one. Even if the ojp ratio is reproducible by pseudo-time- 
dependent models of the type considered here, we must realize 
that the fractional abundance of the cyclic species tells a differ¬ 
ent story. 
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Table 2. Dominant reactions, branching ratios, and rate coefficients for the C 3 H 2 ortho/para problem 


Standard Models 


Alternative Models 



Reactions 

BR-^ 

Default & Sets of Models 1-6 

a"' SM" 

TrialI ^ 

a 

Trial2 

a 

Formation 


C 3 H+ + P-H 2 


F-C 3 H+ 

k\ 

3.30(-13) 

1.0 

• 



0 -C 3 H 2 + c 


0 -C 3 H 2 -H H 

l.Qki 

3.15(-7) 

0.5 

• 




-> 

P-C 3 H 2 -H H 

0.0 k2 

0 

0 






C 3 H H- . 

f 

3.15(-7) 

0.5 




/ 7 -C 3 H 3 + e 


O-C 3 H 2 -H H 

0.5 k2 

1.58(-7) 

0.5 

• 





P-C 3 H 2 + H 

0.5 k2 

1.58(-7) 

0.5 

• 





C 3 H H- . 

/ 

3.15(-7) 

0.5 




O-C 3 H 2 -H H+ 

-> 

0 -C 3 H+ -H H 2 

o.eiki 

5.16(-9) 

0.5 

• 

5.16(-9) 

1/5 x5.16(-9) 



/ 7 -C 3 H+ + H 2 

0.33 k 3 

2.54(-9) 

0.5 

• 

2.54(-9) 

1/5 X 2.54(-9) 

0 -C 3 H 2 + H 3 O+ 


0 -C 3 H+ -b H 2 O 

0.67 k' 

2.48(-9) 

0.5 


2.48(-9) 

1/5 X 2.48(-9) 



/ 7 -C 3 H+ + H 2 O 

0.33 k'^ 

1.22(-9) 

0.5 


1.22(-9) 

1/5 X 1.22(-9) 

0 -C 3 H 2 -H HCO+ 


0 -C 3 H+ -b CO 

0.67 k” 

2.14(-9) 

0.5 


2.14(-9) 

1/5 x2.14(-9) 



/ 7 -C 3 H+ -b CO 

0.33 k” 

1.06(-9) 

0.5 


1.06(-9) 

1/5 X 1.06(-9) 

Destruction 









/ 7 -C 3 H 2 + H+ 


0 -C 3 H+ -b H 2 

0.0 ki 

0 

0 

• 

0 

0 



/ 7 -C 3 H+ + H 2 

1 . 0^3 

7.70(-9) 

0.5 

• 

7.70(-9) 

1/5 X 7.70(-9) 

/ 7 -C 3 H 2 H- H 3 O+ 

-> 

0 -C 3 H+ -b H 2 O 

0.0 k'^ 

0 

0 


0 

0 



/ 7 -C 3 H+ -b H 2 O 

1.0 k' 

3.70(-9) 

0.5 


3.70(-9) 

1/5 X 3.70(-9) 

P-C 3 H 2 H- HCO+ 


0 -C 3 H+ -b CO 

0.0 k” 

0 

0 


0 

0 



P-C 3 H+ -b CO 

1.0 k" 

3.20(-9) 

0.5 


3.20(-9) 

1/5 X 3.20(-9) 

[ 0 , / 7 ]-C 3 H+ h- C 


products 

Jc^ 

1.00(-9) 

0 

• 



[ 0 , / 7 ]-C 3 H+ h- 0 


products 

k' 

4.50(-ll) 

0 




[ 0 , / 7 ]-C 3 H+ h- S 


products 

k" 

A 4 

1.00(-9) 

0 




[ 0 , / 7 ]-C 3 H+ h- Si 


products 

k'" 

K-4 

2.00(-10) 

0 




[ 0 , / 7 ]-C 3 H 2 + C+ 


products 

ks 

4.20(-9) 

0.5 




[ 0 , / 7 ]-C 3 H 2 + S + 


products 

k' 

3.16(-9) 

0.5 




[ 0 , / 7 ]-C 3 H 2 + Si+ 


products 

^5 

3.28(-9) 

0.5 




Interconversion by 








0 -C 3 H 2 + H+ 


P-C 3 H 2 + H+ 

ke 


0.5 


1/4 X 1.30(-8) 

l/4x 1.30(-8) 

P-C 3 H 2 + H+ 


0 -C 3 H 2 -b H+ 

k-6 


0.5 


3/4 X 1.30(-8) 

3/4x 1.30(-8) 

Interconversion by HX^ 








O-C 3 H 2 -H H+ 


P-C 3 H 2 + H+ 

h 


0.5 


1/4 X 7.70(-9) 

1/4 X 7.70(-9) 

P-C 3 H 2 + H+ 


0 -C 3 H 2 -b H+ 

k-i 


0.5 


3/4 X 7.70(-9) 

3/4 X 7.70(-9) 

O-C 3 H 2 -H H 3 O+ 


P-C 3 H 2 -b H 3 O+ 

k; 


0.5 


1/4 X 3.70(-9) 

1/4 X 3.70(-9) 

P-C 3 H 2 H- H 3 O+ 


0 -C 3 H 2 -b H 3 O+ 

k'_j 


0.5 


3/4 X 3.70(-9) 

3/4 X 3.70(-9) 

0 -C 3 H 2 -H HCO+ 


P-C 3 H 2 -b HCO+ 

k" 


0.5 


1/4 X 3.20(-9) 

1/4 X 3.20(-9) 

P-C 3 H 2 -b HCO+ 


0 -C 3 H 2 -b HCO+ 

k-2 


0.5 


3/4 X 3.20(-9) 

3/4 X 3.20(-9) 


Note: a{-b) - ay. 

“ The interconversion rates are obtained as discussed in Section l2^ 

* Competitive protonation rates are suppressed by factors of 5. 

Branching ratio between ortho and para modifications obtained by angular momentum algebra ilokal2nn4 . 
k- a (T/300)“^ exp (-y/T); a is in units of cm^ s“'. 

SM (Simple Model): the reactions marked with • are used in the model of iMFKl: with the generic forms HX^ and Y. See 
Section 3.3. 

In our standard approach, the branching fraction to form C 3 H is set equal to that for C 3 H 2 except for models where it is a 
free parameter. 













